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Abstract 

oo 

We study the Kapteyn series ^^t"J„ (nz). We find a series repre- 

n=l 

sentation in powers of z and analyze its radius of convergence. 
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1 Introduction 

Series of the form 

oo 

[{n + u)z], z/ G C (1) 

n=l 

where J„ (■) is the Bessel function of the first kind, are called Kapteyn series. 
The series is convergent and represents an analytic function throughout 
the domain in which 17.3] 

n{z)< lim |«„r^"+"\ (2) 



n^oo 
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where 



;exp [a/1 — 



1 + 

Kapteyn series first appeared in Lagrange's solution [7] 



(3) 



oo 

E = M + 'i^-Jn {en) sin [Mn) 



(4) 



ra=l 



of Kepler's equation j3] 

M = E - esin(M). 

The solution was independently discovered by Friedrich Bessel in 
where he introduced the functions which now bear his name. 

Kapteyn series were systematically studied by Willem Kapteyn in his 
article 0, where he proved the following expansion theorem. 

Theorem 1 Let f{z) be a function which is analytic throughout the region 

Da = {zeC \ n{z) < a} , 

with a < 1. Then, 

oo 

Q^nJn {nz) , Z G Da 



n=l 



where 



ttn = ^ y ©n (z) f{z)dz 



and the path of integration is the curve on which VL[z) = a. The function 
B„ (z) is the Kapteyn polynomial defined by 



ejz) 



©n (Z) 



I - I 

lk^{n- 2kf {n-k- 1)! fnz\^k-n 



E 

k=0 



2 J 



, n > 1, 



where [-J denotes the integer-part function. 
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Proof. See 17.4]. ■ 

As a corollary |H], one finds that if the Taylor series for f{z) is 



then, 



oo 

,.-2' 

n=0 



Itl 9 

1 ^ (n-2A;)^(n- A; - 1)! 

ZZ^ ,, ,„xn-2fc+i «""2fc, ri>l. (5) 



an ^ 

fc=o ^' (i) 



In IHj, Kapteyn established the fundamental formula 

^ oo 

— = 1 + 2^J„ (nz) , Q{z) < 1, (6) 



1 - z 

n=l 



using the integral representation |H] 



27rz J \ s I "5 

where the Hankel contour Ti. encircles the origin once counterclockwise. 
The purpose of this paper is to generalize (0) by means of the function 

oo 

F{z, t) = (nz) , zeC, teR, (7) 

n=l 

so that (jni) corresponds to the particular case t = l.In Section 2 we write 
F{z,t) as a series in powers of z 

oo 

F{z,t) = Y,Mt>'' (8) 

n=l 

and find the coefficients An{t). Although one could replace = in (jSj) 
and solve the problem 



we take a different approach that only uses the differential equation satisfied 
by the Bessel functions. In Section 3 we analyze the radius of convergence 
of (jHl) for different ranges of t. Finally, in Section 4 we solve the inversion 
problem (0) for arbitrary an- 

2 Power series 

Since the Bessel function Jn{z) is a solution of [Hj 

z'^J'^ + zJ'^ + (2^ - n^) Jn = 0, 
the function yn{z) = Jninz) satisfies the ODE 

z'y'n + zy'n + {z^ -l)yn = 0. (9) 

We also observe that 

f{r)" + t{r)' = nH''. (10) 

Using © and ^ in ((Tj), we have 

z^F,, + zF, = (1 - z^) {t^Fu + tFt) , (11) 

where the subscripts denote partial derivatives. 
Replacing (jH)) in (fTT|). we get the equation 

t'A: + tA'^-n^An = t^A:_^ + tA'^_2, n>l, (12) 

where we define An = for n = —1, 0. 
From and (jH)), we have 

00 

= F(;2,0) = ^A„(0)^" 

n=l 

and, from we see that 

00 _ 00 

n=l ^ ' n=l 

which together imply that 

AniO)=0, A„(l) = ^, n>l. (13) 



Solving (fT^ with the boundary conditions (fT^. we obtain 
t 9 1 

Ai(t) = -, A2(t) = — , Agft) = —t^ 1 

L\ J 2 ^ ^ 2 ^ ^ 16 16 

. / X 2 4 1 2 . . X 625 , 81 o 1 

AJt) = -t^ - -^^ ^5 m = + 1, .... 

^ ^ 3 6 ^ 768 256 384 ' 

Thus, the function An(t) is a polynomial in t of degree n. Writing 

n 



k=0 



and using (fT^ . we get the recurrence 

(n^ - A;2) CI? = -A;=^C^-^ for A; = 0, 1, . . . , n - 2 



e) CI! = 0, for A; = n - 1, n. 



[n 



To solve (fT3j) . we recall that the double factorial function satisfies 

(n - k)\\ = {n-k){n-k- 2)!!, {n + A;)!! = (n + A;) (n + A; - 2) 
which gives 

/ 2 _ l2\ \ _ \ 

' {n-k)\\ X {n + k)\\ ~ {n~k-2)\\ x (n + A;-2)!!' 
This suggests the solution 

cos (^vr) A:" 

Cfc = 7 ,,,, ' / TTTv 0<k<n. 

" {n- k)\\ X (n + A;)!! ~ ~ 

Therefore, 



" cos (^tt) A:" 



or, after rearranging terms, 

LfJ 



n 



It is clear from (fT7|) that A„(0) = 0. To show that = |, we use the 

identity H 0.154, #5] 

E(-l)'Q(A^ + «r = (-irn!, (19) 

fc=0 ^ ' 

where a is arbitrary. From ()18|) we have 
Setting t = 1 in (j^ . we have 

fc=0 ^ ^ 

and from (|T9|l . we conclude that 2y4„(l) = 1. 



3 Radius of convergence 

We shall now find the radius of convergence Rif) of the power series (jH)). From 
((7j) we observe that — > as t — cxd and R{f) — oo as t — * 0, while (jH)) 
gives -R(l) = 1 (see FigureC}. Since from (fT8|) we have A„(— t) = (—1)" Anif)^ 
we shall limit our analysis to t > 0. 

When \, the largest term in the sum (fTH|l corresponds to A; = 0, and 
therefore 



1 (nV 



^n{t)---A-\ , t^OO. (21) 



Using Stirling's formula PP 



we obtain 

'''' 

Thus, 

1 26"-*- 

R{t) = lim \An{t)\'" , t oo. (23) 

n— >oo t 
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Figure 2: A sketch of In [|A„(0.1)|]. 



When f < 1 , (jl8|) is highly oscillatory and there is no approximation like (1211) 
valid in this range. 

To find a formula for R{t) valid for all t, we observe that, for a fixed value 
t > 0, we have 

ln[|yl„(t)|]^/?(fK (24) 

for some function (see Figure |21). 

Thus, we consider an asymptotic expansion of the form 

An{t)^^9nit)[ijit)f, n^oo, (25) 

where On{t) denotes the sign of v4„(t) and ipi^l) = 1. Replacing in (fT^ 
and setting 

e'„it) = o:it) = o, 9'„_,it) = e:_,it) = o, 
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we obtain, to leading order in n, 

{t^'f _ _ Q^^^^A ^ Q_ (26) 

Since 

dn = ±1, dn-2 = ±1, 

and if){t) is increasing for t > 0, we get 

^' = (27) 
Solving (f?rj) subject to '0(1) = 1, we obtain the implicit solutions 

ln(f) + M!i^I±ill ^^2(^) + i_arcsinh [^{t)] - v^+ln f + l) = 0, 

(28) 

for < t < 1 and 

3 

in(t) - ^"^'^^^'^ + mvnt) - 1 - In [m + vnt) - 1] = o (29) 

for t > 1. 

Although we cannot solve (^Hj) and (j^^ exactly, we can consider the 
limiting cases as t — and t ^ oo. Since 'ilj{t) — > as t — 0, we obtain from 

(121 



1 



t ^ 0. 



-ln(t) + y2 + ln(V2- 1) 
When t — > oo, we have ip{t) — > oo, and we get from ()29j) 

, , , e 

which agrees with (j22j) . 

Since R(t) = we have from (^Hj) and (j^ . after exponentiating 



(30) 



(31) 



i?(t)exp Vl + i?2(t) 

e-^ (1 + V2 ) ^^=^^ 

V / 1 + v/i + i?=^(t) 



■t = l, < t < 1, 



(32) 
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and 



R{t) exp [v^l - R^{t) 



t=l, t>l. 



, . ^. (33) 

1 + Vi-^'W 

In Figure El we graph the solutions of ip^. (fH^ and the approximate value 

of Rit) given by 1^500(01"^- 

Using with a„ = t", we conclude that the Kapteyn series ((Zj) converges 
for those t G M and 2; G C such that 



;exp [vl 



< 1. 



(34) 



Replacing z = rd^ in ()34|1 . we find that the minimum value of r corresponds 
to = ±|. Thus, for t > 0, the Kapteyn series ((7|) will converge inside the 
circle \z\ < r{t), with 



r{t) exp a/1 + r(t)2 



1 + v/l+r(t)2 



(35) 



As it was observed in [S] , the circles of convergence of the Kapteyn series ((Tj) 
and the power series (jH)) are not equal, the former being slightly smaller than 
the latter. However, the radius r(t) given by ()35|) and R{t) are asymptotically 
equal as t ^ and t ^ 00 (see Figure Ej). 

We summarize our results in the following theorem. 

Theorem 2 The Kapteyn series 

F{z,t)=J2t^Jrr 



nz) 



n=l 



converges for those t G M and z E C such that 

;exp [Vl — -2 



< 1. 



1 + 

The function F{z, t) admits the power series representation 

00 

F{z,t) = Y,A4t)z-, 

n=l 
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Figure 4: A comparison of r{t) (solid curve) and R{t) (ooo). 
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with coefficients 



An{t) 



-1) 



k=Q ^ ^ 



2/ 



and radius of convergence R{t) defined by the implicit equations 

-t = l, < |t| < 1, 



exp ^/T+mt) 

6-^(1 + V2) )= 



■t = l, |t| > 1. 



1 + v^r^^ 

4 The general problem 

We shall now consider the problem of finding a„ in terms of a„, where 

oo oo 



n=l 



n=l 



for arbitrary a„ and a„. Rewriting the power series of J„ (nz) [S] 

(—1)-' /n2;\"+2i 



j=o 



+ V 2 J 



in the form 



we have, 



oo 



n=l 



n=l 

oo k 

EE 

fc=l n=l 



a„ 1 - 



n\«cos [f (A;-n)] (f) 



k—n 



(^)!(^)! 
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Therefore, we obtain from (j36p that 



k 

COS 

a. 



n=l 



or 



^ cos [| {k — n)] 



^ {k-n)\\{k + n)\\ ^ ' 



n=l 



In particular, setting a„ = in (jSTj), we recover (fT7|) . 
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